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MODELING OF THIN-FIUl GaAs GROWTH 
By 

John H. Helnbockel'*' 

SUMMARY 

A solid-on-solid Monte Carlo model is constructed for the simulation of 
crystal growth. The model assumes thermally accommodated adatoms Impinge 
upon the surface during a At time interval. The surface adatoms are 
assigned a random energy from a Boltrmann distribution, and this energy 
determines whether the adatoms evaporate, migrate, or remain stationary 
during the At time interval. For each addition or migration of an adatom, 
potential wells are adjusted to reflect the absorption, migration, or 
desorption potential changes. 

INTRODUCTION AND STATEMENT OF PROBLEM 

Numerous methods have been applied to obtaining thin film, single 
crystals of GaAs, including free-standing wafers, peal films removed from a 
single crystal substrate, and films grown on light-weight substrates. One 
of the most promising methods is a version of the last technique called 
"graphoepltaxy . It is generally known that overlayers of crystalline 
materials deposited upon smooth microcrystalline substrates tend to be more 
or less randomly polycrystalline. The absence of long-range order in the 
microcrystalline substrate is reflected in the absence of long-range order 
in the overlayer. The basic concept of graphbepita^ is that, by intro- 
ducing an artificial surface relief structure having a long-range order on a 
microcrystalline substrate, long-range order can be Induced in an overlayer. 
In other words, a crystalline film can be grown on a microcrystalline 
substrate. 


^Professor, Department of Mathematical Sciences, Old Dominion University, 
Norfolk, VA 23508. 


In order to be able to study the dynamics of crystal growth on grat- 
ings > a more detailed model of crystal growth Is needed. The following Is 
the development of a solld-on-solid (SOS) model for the simulation of molec- 
ular exchange between the solid and vapor phases whereby dynamics are intro- 
duced into the model by a random deposition, migration, and desorption of 
adatoms at random surface positions. 
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lattice constant (nm) 
size of square array 
Integer position 
height at position I 

potential energy change at site I caused by adatom deposition (eV) 

(1 ■ 1,...,8) potential energy changes at sites neighboring I (eV) 
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pocentlal at site I (substrate) (eV) 

Boltzmann constant (8.6 x 10**®) (eV/*K) 

absolute tei^perature C*K) 

scale factors for Mie potentials 

U .(ijk) potential energy on (ijk) face for cluster of n^adatons (eV) 
nCuUSTER 

B .. median energy for Boltzmann distribution 

median 

Integers associated with (ijk) orientation 
AUjj^Ujj-Uo surface diffusion activation energy (eV) 

AU^»0-Uq evaporation energy (eV) 

reference temperature (®K) 

R uniform random number 


K 
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frequency factors 
diffusion coefficient 
mean adatom displacement 
deposition rate 
evaporation rate 
vertical growth rate 
supersaturated vapor pressure 
equilibrium pressure 


MODELING OF GaAs CRYSTAL GROWTH 


Structure of GaAs 

By definition, a crystal is any solid with an ordered arrangement of 
its atoms. The cyrstal structure of GaAs can be thought of as a lattice 
with repeated units of face-centered cubic (fee) cells. Consider two types 
of fee cells, one with Ga atoms and another with As atoms. Imagine these 
two fee cells coalesced with the same origin at one corner. Now slide the 
origin of the As cell along the vector (a /4)[lllJ where a Is the lattice 
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constant. The resulting fee arsenle cell has its origin displaced to the 
point ( a^/4,a^/4,a^/4) with respect to the gallium fee origin. These inter- 
penetrating fee cells form a pattern which constitutes Che GaAs lattice 
structure. The value of the lattice constant is a^ ■ 0.56417 nm. 

SOS Model 

The solld-on-solld (SOS) model is a square array of "blocks” stacked to 
form columns. This type of model can be represented by a square array of 
integers, where each Integer represents the number of "blocks" in each col- 
umn perpendicular to some reference plane. The "blocks" can represent atoms 
or molecules which are being stacked. Hereafter "blocks" will be referred 
to as adatoms. The SOS model is a generalization of the terrace -ledge-kink 
model (ref. 1) , and Che surface elements are those adatoms at the tops of 
their columns (see fig. 1). It is assumed that adatoms can be added (or 
removed) only to (from) positions over other solid adatoms. 

In the SOS model the columns are built up upon an M x M square array by 
randomly placing adatoms upon the array and allowing these randomly depos- 
ited adatoms to either condense, evaporate, or migrate. The M x M array can 
be replaced by a single subscripted array so chat only one random number is 
needed to select an arbitrary column (ref. 2). This one-to-one correspon- 
dence is Illustrated in table 1. In the table, note that additional rows 
and columns have been added to the array. These augmented rows and columns 
will be used to eliminate edge effects and will be assigned periodic bound- 
ary conditions. We let H(I) denote the height of a column at some arbi- 
trary position I where n + The periodic boundary condi- 

tions imposed upon the heights are Chat they satisfy the conditions: 

(1) Row 0 is Che same as row M and 

H(I) - H(I + m2 ), I - 10 n + 1; 

(2) Row 1 is Che same as row M + 1 and 

H(I) - H(I + m2), I + n,..., n + M -»• 1; and 

(3) Column 0 is the same as column M and column 1 is Che same as col- 
umn M -f 1 when the columns are displaced by one unit. 
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ADATOM ^ / STEP OR jOG 



ADATOMS 


Figure 1. SOS model (K ■ kink site). 


Table 1. Integers associated with row— column positions of M 
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Thus» to each position I (n + 1 _< I _< n + ), we can assign a height 

H(I) which represents the number of stacked adatoms. Also associated with 
a general position I there are neighboring positions which have the posi* 
tlon indices given in table 2. The SOS model can be described as an array 
of interacting columns of Varying Integer heights. The surface adatoms are 
at the tops of columns. 

The term "epltai^'' means "an arrangement on" and is used to denote the 
growth of one substance upon the crystal surface of a foreign substance. 

The term "auto-epitaxy” is the oriented growth of a substance onto Itself, 
and "heteroepitaxy" is the growth over another material. Obviously, hetero- 
epitaxy becomes autoepitaxy after one layer of growth has been deposited. 

We will use the SOS model to simulate epitaxial growth. 

Potential Energy of Adatoms 

The "rules" by which tb© columns of the SOS model Interact will be 
governed by the following ideas relating to potential energy and the poten- 
tial energy changes which are associated with the adsorption, migration or 
desorption of adatoms. Vfe define the following energies: 


Uq ■ Uq(I) potential energy at site I due to surface bonding 

and crystal structure 

potential energy change at site I due to the adsorp- 
tion or desorption of adatoms (the same for all 
sites) 

“U)j^(i»l, . . . ,8) potential energy changes at neighboring sites when 

an adatom is deposited at I (the same for all sites) 


E(I) 


random surface energy associated with site I and 
time interval At 


U(I) ■ U (I) total energy during time interval At. 
+ E(I) ® 


evaporation potential 
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Table 2. Indices associated with neighbors sur- 
rounding a general position I within 
the M >« M array of the SOS model. 

' - — -■ 

I-M-1 I-M I-M.+l 

I - 1 I I + 1 

I + M-1 I + M I + M+1 






migration potential 


All of the above energies are measured in electron volts* 

We will develop a Monte Carlo coo^uter simulation of crystal growth 
(refs* 3-8) by developing "rules" which determine the SOS kinetics of con- 
densation, evaporation, or surface migration of adatoms. Consider first the 
adsorption of an adatom onto the surface at some general site 1* In addi- 
tion' to the change in the potential at I, potential energy changes simultan- 
eously occur at the sites neighboring the central site I. The potential 
energy changes are depicted by the mnemonic "mask" of table 3 which caa be 
placed over the central site I to illustrate what changes must be made in 
the surrounding potential wells. We consider next the desorption of an ad- 
atom from the surface* In this case the potential energies at I and its . 
neighboring sites again change and we use the mask of table 3 with opposite 
signs on the potential changes. 

The case of surface migration is treated as a desorption from a site I 
follotred by an adsorption at a nearest neighbor site to I together with the 
correct potential mask changes associated with each process. A random walk 
is performed to determine nearest neighbor migration sites. 

The Monte Carlo simulation of crystal growth involves a random deposi- 
tion of thermally accommodated surface adatoms during a time interval At. 
These deposited adatoms will change the potential energies at the surface 
sites under consideration. The values assigned to 4*^ -m^,i«l, . . . ,8 

will dictate the changes in potential energies. We assume that the surface 
adatoms have a known a priori surface energy distribution f(E) from which 
we can assign random energies to the surface adatoms* We let 

U(I) - U^(i) + E(I) (I) 

denote the total energy assigned to a surface adatom during the time inter- 
val At. This total energy is the sum of the potential energy 
the structure of the lattice and the random energy E(l) from the assumed 
surface energy distribution f(E). If U(I) is less than some material- 
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Table 3. 

Potential 
with each 

energy adjustments 
site I. 

assi r.iated 




“<Dj 

-0)g 


♦o 


-«5 


-<i*4 

““3 



dependent migration level U , then the adatom w±ll remain at the site I. 

m 

If U(T) Is greater than some materlal*'dependent evaporation level then 

the adatom will evaporate from the surface. For U < U < U , then surface 

m e 

migration Is allowed to occur. During each time Interval At there Is a 
random deposition of adatoms upon the surface. Each of these new adatoms to- 
gether with the other surface adatoms Is assigned a random energy and tested 
to deteirmine If it remains on the surface, migrates, or evaporates from the 
surface. 

This type of model can simulate the atomic or molecular exchange be- 
tween vapor and solid phases of crystal growth. We assume that an adatom Is 
a two-component species of Ga and As. Also, in the Monte Carlo simulation. 

It is assumed that the rate of Impingement of adatoms is independent of the 
neighboring surface configurations. The number of thermally accommodated 
adatoms which arrive on the surface is given by 


P 

(2TrmKT)^ 


( 2 ) 


where P Is the vapor pressure, M Is the vapor adatom mass, K Is 
Boltzmann’s constant, and T Is the absolute temperature. 

The rates associated with evaporation and migration of adatoms depend 
upon the neighboring potentials which must be broken during these processes 
and also upon the random surface energy distribution which Is assumed. The 
fate of a surface adatom depends upon Its total energy U at any Instant of 
time as well as upon where this energy lies with respect to the condensation 
and evaporation energy levels and U^. 

In assuming values to the potential changes ij) and o). ,1«1 , . . . ,8, we 

o 1 

must take Into account the type of crystal structure and orientation we are 
trying to simulate with our SOS model. Consider figure 2 which Illustrates 
the GaAs fee structure. For growth on the (100) face we can set up a corres 
pondence between a central site, the nearest neighbor sites, second nearest 
neighbor sites, and the adatom potential changes for the mask in table 3 

(l.e. , o»i ■ <j»2 » W2 " etc.). Similarly, we can set up the correspondences 
illustrated in figures 2 (b and c) for the (111) and (110) orientations. 
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In figute 2 we oust choose ♦o *<j>2 »^3 such a way that, when the 
first level of adatoms covers the surface, then the potential distribution 
must return to its original value* We will require that adjustments be made 
in the potential energy changes during the transition from heteroepita]^ to 
autoepitaxy. Here we let a negative sign denote an attractive potential* 

By simply adding adatoms to a surface it is readily verified that the poten- 
tial changes must adhere to the rules given in table 4 if after one layer 
the potential energy returns to its initial value* 

We can assign arbitrary values to the neighbor potential changes 

assume a (m^,n^) Mle potential, with < n^, (ref, 9); 

^ 

which is Illustrated in figure 3, Here R* Is the distance at which $ 
obtains Its minimum value of 
various R*/R values* 

The values assigned to the mask potential changes are not necessarily 
the same for the different orientations: for example, the <^i*^2>^3 values 

for each case In table 4 could have different values* 

If we assume a (3,9) Mle potential and an fee crystal, we can calculate 
the potential changes caused by adding adatoms to neighboring sites* If we 
leave the central position empty and fill nearest neighbor sites first, then 
second nearest neighbor sites, followed by adatoms to third nearest neighbor 
sites, and then fill In the central position, we will obtain the curves In 
figure 4* In figure 4 the lower curves are when the site I is Initially 
empty and the upper curves are for the site I Initially occupied by an ad- 
atom* An examination of figure 4 shows that, If we add a random energy to 
the curves Illustrated, then the upper curves represent the potential ener- 
gies associated with groups of adatoms clustering about I* If we consider a 
cluster of adatoms which is smaller than some critical cluster size, then 
this cluster will be unstable in the sense that an adatom at site I will 




Table 5 gives values of 4>/(^ for 
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Figure 2. GaAs cell with growth patterns and potential changes for (100), 
(111), and (110) faces. 
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Table 4. Potential changes for addition of an adatom to an 
arbitrary site. 
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Table 5. Mle potentials 
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CHANGE IN POTENTIAL ENERGY AT SITE I 



Figure. 4. Change in potential energy at site I due to deposition of 

adatoms to neighboring sites (nearest neighbors filled first 
with assumed (3»9) Mle potential). 




have a probability of greater than 1/2 of leaving the cluster. In the case 
where the cluster size Is greater than the critical cluster size, then we 
desire such a cluster to be stable In the sense that an adatom at site I 
will have a probability of less than one-half of leaving the cluster. 


For a (m^,n^) Mle potential, we can express the second and third near- 
est neighbor potentials ^2*^3 terms of the nearest neighbor potential 
, and we can write 


♦2 “ 


(4) 


4'3 “ “ 3^1 

where 012*03 are appropriate scale factors (see, for example, table 5). 
Using the potential masks of table 4, we can calculate the potential ener- 
gies at an arbitrary site I for different crystal orientations and various 
cluster sizes about the central site. These various potentials are given in 
table 6 for arbitrary ,03 ,<J>i values. 

In our computer model we will assign a value to <j>^, the potential 
change due to the addition of an adatom at site I. Then from the (®q» u^) 
Mle potential and the relations • 4j>j^ + 4^2 " 4 qi2 )(^1 (100 case), <})^ 

• + 2^2 ■ (6 + 2o 2 )<^1 (Hi case) and ■ 2(j)i + 2^2 "*■ ^3 * (2 + 2a2 

+ ^3>^l (110 case), we can determine the nearest neighbor potential change 
and from equation (4) the other potential changes for nonnearest neigh- 
bors. In our model we will use a (4,10) Mle potential (refs. 9, 10). 

The data presented in figure 4 and table 5 can be used to make a rough 
limited estimate of the nearest neighbor potential (^x in the following 
sense: we assume that the surface thermodynamics can be attributed to the 

property of clusters. Using the Gibb's free energy of formation of a spher- 
ical cluster of radius r as 


AG - 4ii r^o - A V r^ — In — 

vr Tt 


where r “ radius of cluster, o ■ surface vapor interface energy, V ■ molecu- 
lar volume, T ■ absolute temperature, P ■ supersaturated vapor pressure, 

■ equilibrium pressure (P/P ■ supersaturation). (See figure 5 for a dia- 

e 

gram of potential energy.) Here AG is a maximum at 
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* 2aV 

r ■ r ■ 

KT In |- 
e 


This is the radius of the critical cluster. If a cluster has a radius r 

A 

which is smaller than r , then the cluster will be unstable and an adatom 
will have a probability of greater than one-half of leaving the cluster. 

A 

Alternatively, if the radius of the cluster is larger than r , then an ad- 
atom will have a probability of less than one-half of leaving the cluster. 

These ideas can be Incorporated into our model in a limited sense as 

follows: We let » o denote the evaporation level; and AUg then 

represents the desorption energy of a surface adatom and A^ ■ or 

U « “AU , for AU , the activation energy for migration, we have U - U » 
o e’ m ^ Of m o - 

AU . By using the results of table 6, other potential energies can be 
m 

visualized inside the potential well. Usually, AU ,AU are measured at 

e m 

reference T^, and it will be assumed that these values are temperature 

independent. We also assume that the thermally accommodated adatoms have a 

Boltzmann distribution which has the median value E * KT In 2, where 

median r 

K is Boltzmann's constant. If we assume a critical cluster size and choose 
such that this cluster has an equal probability of growing larger or 
smaller, we have for the different crystal orientations the requirement 
that, when a median energy from the Boltzmann distribution is added to the 
energy of an adatom associated with a critical cluster, we obtain 


U + g 4). + KT In 2 - U 

o ®n r m 

energy of critical median migration 

cluster adatom energy energy 


Equation (5) enables us to determine (j>i such that 50 percent of the 
adatoms will remain at the site I, while 50 percent of the adatoms will 
desorb from the site I. An example which will illustrate the use of equa- 
tion (5) is given in the next section. Also, the limitations of using equa- 
tion (5) to assign a value to (pi will be discussed. 
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Table 6. Potential energies at arbitrary site I. 


Cluster 

Size n 100 


1 Uq + (4 + 402 )*i 

2 + (3 + 4o2)'|>i 

3 Uq ♦ (2 + 4a2>^i 

4 1^0 * 4a2>^i 

5 Uq + 402^1 

6 Uq + 302^1 

7 Uq + 202^1 

8 + Q«2«^1 


Potential at Site I 

gnh 

m . 


+ (6 + 202)^1 
Uq + (5 + 203)^1 
+ (4 + 2 o 2 ) 4 »i 
Uq + (3 + 2 o 2 ) 4 »i 
Uo + (2 + 2a2)4'i 

Uq + (1 + 2o2)^i 

Uq + 202 ♦! 

Uq + «2^1 


110 


Vq + (2 + 2 o 2 + 403)411 
Uq (1 + 202 + 403)411 
Uq + (202 + 403)411 
Uq + (c«2 + 4 o3>^i 
Vo + 4034)1 

Uo + ^3h 
Uo + 2034*1 
Uo + “a*! 
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Energy Distribution and Random Walk 


In our model of crystal growth we assume that the thermally accommo*' 
dated adatoms possess a Boltznunn distribution for their surface energies* 
This distribution can be i^ixpressed as 



It has the cumulative distribution function: 


F(E) - / f(E)dE - 1 - exp (7) 

We can generate a random variate E from this Boltzmann distribution by 
using the Inverse function associated with the cumulative distribution [eq. 
(7)]. For R, a uniform random number between 0 and 1, the random energy Is 
given by 

E - -KT In (I “ R) (8) 


The mean of this distribution Is KT and the median Is given by 

E - KT In 2 (9) 

median 

We will choose the nearest neighbor potential by assuming a criti- 

cal cluster size (refs* 11, 12) from table 6* For example, if we choose 3 
as the critical cluster size, then we desire that an adatom added to site I 
will have a probability of less than one-half of leaving the cluster* For a 
(100) orientation, we require that 

U (100) + E - U 

3CLUSTER median m 
or 

Uq + (2 + + KTj. In 2 - 
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to satisfy 


This requires that we choose 


- U - KTr In 2 
m o 

(2 + 4o2 


AU - KTr In 2 
m 

' (2 + 402 ) 


( 10 ) 


where Is the surface diffusion activation energy. For this choice 

of we will then have the following Inequalities: 

for cluster sizes greater than three, 

U (100) + E < U , n > 3 

nCLUSTER median m 


for cluster sizes less than three, 


U (100) + E > U , n < 3 

nCLUSTER median m 

The appropriate choice of then produces the stability and instability 

characteristics associated wlt.h various size clusters. 

In general, we can replace the denoniinator In equation (10) by the more- 
general e:q>ression g^ • (n^+ mjoi + liag) where ni,m^,ti are Integers asso- 
ciated with different crystal orientations. This produces the more general 
expression 

AU - KT In 2 

( 11 ) 

For the nearest neighbor potential c age associated with a critical 

cluster and surface orientation of table 6. In equation (11) the Mie poten- 
tial Is assumed, and second, third nearest neighbors' potential changes are 
given by 

♦2 » a2<fri , <^3 “ 03^ 1 
and T^ is a reference temperature. 

Ue can thus summarize the Monte Carlo simulation model as follows: 

(1) We assign an Initial geometry and characteristics to a sub- 
strate. In addition to geometry considerations, we assign a 


uniform pocential (U ) together with evaporation (U } and ml' 

08 es 

gratlon (U ) levels, 
ms 

(2) HeteroepitajQT Is assumed to occur during deposition of the first 

layer. A flat uniform surface would have potential changes from 

a uniform value of U everywhere on the substrate to a value 

08 

of U everywhere after the first layer covers the substrate. Be- 
o 

cause of this change from heteroplta^ to autoepitaxy , the Ini- 
tial layer potential mask Is adjusted by a value of 
from which we can obtain the heteroepitaxy-to-autoepltaxy poten- 
tial change. 

(3) We assume Mle potentials for both the substrate and growing mate- 
rial. The nearest neighbor potential change is determined 

by AU^ ■ U - U , which Is the surface diffusion activation 
D m o 

energy., temperature, critical cluster size, and surface orien- 
tation. The second and third nearest neighbor potential changes 
determined by the Mle potential scale factors oc 2 ta 3 > 
and <^2 “ > ♦s “ “3^1 then construct the appropriate 

potential masks depicted in table 4 for both the substrate and 
the growing materials. 

(4) During a time interval At we let adatoms impinge upon the sur- 
face and update the potentials U^. 

(5) During the At interval, each of the surface adatoms is assigned 

a random energy (E) from a Boltzmann distribution and the total 
energy at each adatom surface site is tested. For U ■ + E, 

the total energy, we let the adatom: 

evaporate if ^ ^ 

migrate if < U < Ug; and 

remain stationary if U < 

During the heteroepita:w, we use the values U ,U to deter- 
^ " es ms 

mine the fate of an adatom, and for autoepita:qr we use the val- 
ues U and U . 
e m 

(6) We record and update all potential changes and height changes 
which occur during the growth process together with other sta- 
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tlsclcal measures of crystal growth to which we may want to as~ 
sign a quantitative value* 

(7) The process continues until some stopping condition is satis- 
fied. Some stopping conditions are that (a) layers exceed a 
certain value or (b) time exceeds a certain value* 

During the surface migration of an adatom, we allow it to random walk 
to a nearest neighbor site only* Only those nearest neighbor sites which 
are unoccupied and at the same or a lower level are considered first for the 
random walk* If all the nearest neighbor sites on the same level are fil- 
led, then the adatom is allowed to random walk to a nearest neighbor site by 
hopping up to the next highest level* The flow chart of the computer sinu- 
latlon is illustrated in figure 6* 

The limitations in calculating <|)^, rather than assigning a value to 

(which is more general), are as follows: let ■ KT^ In 2, then from 

equation (11) we can construct the following linear relation * m(^U - C ) 

o m o 

where the values of the slope m are given for different cases: 


ciSSUHh* 


100 


Orientation 

111 


110 


2 

m 

■ 4(1 + 02) 

m 

« (6 + 2tt2) m ■ 

(2 

+ 2o2 

+ 403) 



3 + 4 o 2 


5 + 202 

1 

+ 2o2 

+ 4tt3 

3 

m 

- 4(1 + tt2) 

m 

« (6 + 2tt2) m ■ 

(2 

+ 2o(2 

+ 403) 



2 + 4 o 2 


4 + 2tt2 


2o2 + 

403 

4 

m 

4(1 + 02) 

m 

6 + 2 a-y 

m ±. Q ■ 

ii 

+ 2a 2 

+ 403) 


1 + 4 o 2 


3 + 2tt2 


02+ 40^ 


These relations imply that for AU large enough, U + <^ >0; conse- 

m o o 

quently these adatoms will evaporate from the surface* Clearly, we do not 

want A to be dependent upon the choice of AU * As a consequence, the 
o ® 

cluster interpretation still exists, yet we must assign a value to <|)^ which 

is independent of AU * 


( 12 ) 
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Figure 6. Flow chart of model. 
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ADDITIONAL CONSIDERATIONS AND ASSUMPTIONS 


With reference to figure 7, we have the following situation: When 

material is deposited upon an ideally flat substrate, the potential wells 

will change from U to U after one layer of material covers the sub- 

os o 

strate. In order to accomplish this change, we alter the potential mask of 
table 4 by replacing ^ ^o ” ^^os first-layer locations. 

For layers greater than one, we revert back to the masks of table 4. 

The interatomic potential possesses the general characteristics of flg- 

* It 

ure 7. If we assume that « B 2 * interatomic spacing, then 

the first-layer potential adjustment will be as above. If, how- 

ever, there is a misfit dislocation in the atomic spacing of the substrate 

and deposited material, then the first-layer potential adjustment will be in 

* A 

terms of the different potentials at R^ and % and will vary accordingly 
upon the misfit distance. The first-layer potential adjustments will, in 
general, have to be corrected as we move along the misfit direction as sug- 
gested by figure 8. This will require much additional bookkeeping of poten- 
tial values and changes with distance. 

Imagine that we move an adatom along the flat surface. For each dis- 
placement S, measured from some fixed reference, the adatom will experience 
a potential attraction or repulsion which is determined by the Mie poten- 
tial. For each value of S ■ 0,ai , 2ai , . . . , there will be an associated po- 
tential change which can be calculated for each position. This would require 
that the potential masks of the first layer would each have to be adjusted 
based upon the assumed geometry and calculated Mie potential values. 

For a heterogeneous surface, the evaporation potential changes will 

vary from site to site as will the nearest neighbor migration levels. If we 

let U denote the potential well at a position i, j and denote by AU 

ij e 

the average evaporation energy, and AU the average migration energy, then 

m 

associated with each site potential U . there are nearest neighbor migra- 

I » j 

tion levels AU , , ,AU each of the nearest neighbor directions, as 

m 1 m, / 

Illustrated in figure 9. These additional concepts will require a more de- 
tailed computer model to keep track of all the potential changes and updates. - 


27 



Figure 7. Potential energy diagrams for growing material and 
substrate. 
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Figure 8. Misfit dislocation approximation. 
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Figure 9. Heterogeneous surface potentials. 


We begia our study with the oversimplified model flow charted in figure 

6 where we assume a (4,10) Mie potential and let the critical cluster size 

determine and consequently For future studies we will assign 

values to ^ independent of AU values and try to determine the feasi- 

o m 

bility of extending the concepts to an A-B model. The following section 
sunmarizes the processes which the model considers. 


SUMMARY AND DESCRIPTION OF GROWTH PROCESS 


We simulate the number of adatoms impinging upon an area of the 
substrate during a time interval At. These adatoms impinge at random 
positions on the surface and will be absorbed, migrate, or desorb depending 
upon their energy states during the time interval At. The deposited 
adatoms will form clusters (refs. 10, 13). The sizes and shapes of these 
clusters will be determined by the energies of the potential wells, the 
deposition rates, the migration levels, and the surface temperature which 
affect the surface energy Boltzmann distribution. 

When conditions are such that surface diffusion controls the growth, 
the clusters that form will increase in size and eventually coalesce with 
other dusters. After the clusters coalesce, the islands of clusters will 
have empty spaces- between them which will fill in and disappear by the ran- 
dom deposition of adatoms. 

The surface energy is an important mechanism which controls the surface 
diffusion and mass transport of surface adatoms since this diffusion process 
enables clusters to form, then grow and finally merge with other clusters. 

We assume that thermally accommodated adatoms are described by a Boltzmann 
distribut ion 


f(E) 


i- exp (^\ 
KT \KT/ 


E > 0 


(13) 


which has a mean E • KT, mode of median of 2. 
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Note that for the average lifetime of a monomer adatom on the 

surface before It desorbs from the surface, we may write 

T « exp — £ (14) 

s KT ^ ' 

Where AU^ is the evaporation energy from the surface to the vapor and 

is a frequency factor on the order of 10”^^ to 10”^^ seconds at room temper~ 

ature* Letting D denote the monomer diffusion coefficient we write 
8 

5 AU 

Dfi " <*^o ®*P (15) 

where AU is the activation energy necessary for an adatom to jump a lat- 
tice distance d to a neighboring site and u is the frequency factor for 
o o 

the jump* For diffusion-dominated phenomena, the adatoms form clusters by 
way of collisions on the surface* The value of AU^ depends upon where the 
monomer is located (l*e* kink site, edge, free adatom, etc* )* We will 
assume that AU^ is associated with a kink site energy change* From the 
Einstein relation 

Vs • (16) 


where X 
s 


is the mean adatom displacement, we can write 



(17) 


Now if diffusion is the dominant mechanism for growth, then X > d and 

5 O 

we require that AU^ > AUjjj* 

Also associated with the growth mechanism are the deposition rate- and 

migration range and how these quantities interact* For example, with a low 

deposition rate and large values for AU and small AU - there will be a 

e m 

high probability for surface migration of adatoms deposited* This surface 
diffusion is a random walk process which continues until the adatom is 
absorbed by the surface or else desorbs back to the vapor* At the other ex- 
treme, if the deposition rate is much greater than the diffusion rate, then 
the adatoms are being deposited before other adatoms can migrate away by 
surface diffusion* The large deposition rate creates the interaction of 
adatoms and the resultant cluster formation* 
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The size of the crystal formed will depend upon: 

(1) the surface mobility of the adatoms; 

(2) the smoothness of the substrate as a growing surface (i*e* » the 
rougher the surface, the more nucleation sites available and the 
higher the probability of anisotropic column growth. Also, 
there will be a reduced surface mobility of adatoms); 

(3) the deposition rate and the kinetic theory of Interaction and 
chemical bonding of adatoms on the surface which will tend to 
increase the average energy of surface adatoms (l.e. for those 
substrates where there is a strong interaction between adatoms, 
then the adatom surface mobility will be reduced) ; 

(4) the activation energies for surface diffusion and evaporation; 
and 

(5) the orientation of the critical clusters and their density at 
the onset of nucleation. 


If the surface energy distribution Is such that at low temperatures the 
adatom surface migration Is small, then the crystal growth Is mainly from 
condensation on random surface sites. If is the deposition rate and 

the evaporation rate from a surface section which has adatoms and d^ 

Is the depth of a layer, then 


■■o <«d - 


Is the average growth rate. 


(18) 


If the substrate material Is amorphous In structure,- then the potential 
wells will not be periodic in nature. This will decrease the surface mobil- 
ity of adatoms and the impinging adatoms will tend to condense In the neigh- 
borhood of impingement. 


The epitaxial relationship between the substrate material and the con- 
densatlng material, in regards to preferential oriented absorption, depends 
upon the relative potential wells and consequent reaction between the con- 
densate and substrate adatoms. It Is assumed that this? process Is random In 
nature and that the higher index faces grow more rapidly than the lower In- 
dex faces and will eventually disappear as is suggested In figure 10. 
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PRELIMINARY COMPUTER RESULTS 

A computer program following the flow chart in figure 6 was constructed 
Time Intervals of At ■ 0.025 and At <■ 0.01 were utilized for the deposition 
of adatoms, where deposition rates of 200 and 400 adatoms per second were 
used (1 Mm/hr * 0.2778 nm/sec » 200 adatoms/sec for a 20 x 20 array). 

During each At time interval the surface adatoms were assigned a random 
energy from a Soltzmann distribution and the poten'~tlal energies were 
adjusted by using the potential masks of table 4. The potential changes in 
table 4 were approximated by assigning a value to AU and then calculating 
the nearest neighbor potential change from the re-lation of equation 

(10). The other neighboring potential changes were calculated by assuming 
a (m ,n ) Mie potential where <^2 “ 02<^i ^3 " 03(^1* ’Ihese potential 

changes then produced ^ ~ t which depends upon the particular orienta*- 

tlon being modeled* 

Some preliminary results obtained from the computer simulations are 
given in the following figures. In these figures, the initial substrate 
geometry was assumed flat with a ledge of height one along the top row of a 
20 X 20 array. Figure 11 Illustrates the graphic display option of crystal 
growth that is available* Note in the center of this figure the critical 
cluster of size 4 for a (100) orientation of crystal growth. Figures 12(a) 
and (b) illustrate the vertical growth rate vs. time for various orienta- 
tions and migration levels. Note that for AU large (near AU ) there is 

m e 

more evaporation of adatoms from the surface and consequently there is a 
lower growth rate. Note also that the (111) orientation has a larger growth 
rate than the (100) or (110) orientations. 
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Figure 11. Graphic display of crystal growth (100) 
orientation. 
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